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ABSTRACT 

High-resolution observations of the solar photosphere have identified a wide variety of spi¬ 
ralling motions in the solar plasma. These spirals vary in properties, but are observed to 
be abundant at the solar surface. In this work these spirals are studied for their potential as 
magnetohydrodynamic (MHD) wave generation mechanisms. The inter-granular lanes, where 
these spirals are commonly observed, are also regions where the magnetic field strength is 
higher than average. This combination of magnetic field and spiralling plasma is a recipe for 
the generation of Alfven waves and other MHD waves. 

This work employs numerical simulations of a self-similar magnetic flux tube embedded 
in a realistic, gravitationally stratified, solar atmosphere to study the effects of a single mag¬ 
netic flux tube perturbed by a logarithmic velocity spiral driver. The expansion factor of the 
logarithmic spiral driver is varied and multiple simulations are run for a range of values of the 
expansion factor centred around observational constraints. 

The simulations are analysed using ‘flux surfaces’ constructed from the magnetic field 
lines so that the vectors perpendicular, parallel and azimuthal to the local magnetic field vector 
can be calculated. The results of this analysis show that the Alfven wave is the dominant wave 
for lower values of the expansion factor, whereas, for the higher values the parallel component 
is dominant. This transition occurs within the range of the observational constraints, meaning 
that spiral drivers, as observed in the solar photosphere, have the potential to generate a variety 
of MHD wave modes. 
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1 INTRODUCTION 

The solar atmosphere is a highly dynamic, often unpredictable and 
turbulent environment. It also has a direct impact on the Earth, with 
events such as Coronal Mass Ejections causing geomagnetic storms 
which can disrupt systems such as satellites and power grids. The 
outermost layer of the solar atmosphere, the corona, is observed to 
be heated to millions of degrees Kelvin. This hot plasma requires a 
constant energy input to prevent it from cooling, and the mechanism 
by which this energy is transferred into the corona is a subject of 
intense study. The source of the energy is in the photosphere and the 
internal deeper regions of the Sun. 

The solar photosphere is a highly dynamic region of the solar 
atmosphere, with hot plasma rising up from the convection region, 
radiating and sinking back down within inter-granular lanes. Com¬ 
bined with this are multi-scale magnetic fields which intersect the 
photosphere and one outcome is the generation of magnetohydro¬ 
dynamic (MHD) waves by various driving motions and at different 
frequencies. These MHD waves are generated in a variety of vertical 


or near vertical magnetic structures, which connect the different 
layers of the gravitationally stratified solar atmosphere. This yields 
a potential mechanism for energy transport vertically through the 
solar atmosphere, along these magnetic structures, which has been 
widely studied as a potential solution to the coronal heating problem. 
This, MHD wave heating of the solar atmosphere, has been studied 
analy tically (g.g.|Andries et af] 2009, Wang 2 01 1|>, observationally 


(g.g.|Bogdan & Judge 2006, Kobanov et al.|2006||Morton e t al.|2012| 

|Jess et al.|2009||Taroyan & Erdelyi|2009||Dorotovic et al.|2014|> and 


numerically ( e.g . Bogdan et al.|2003||Hasan & van Bal legooijen 


2008 ; |Scullion et al.|2011||Fedun et al.|2011||Vigeesh et al.|2012| 


Wedemeyer-B ohm et al.|2012| > 


This work, as a follow-up to Mumford et al. (2015), investi¬ 
gates the effect of logarithmic spiral-type velocity drivers in the solar 
photosphere and their properties as MHD wave generation mech¬ 
anisms. Mumford et al. ( 2015) studied five different photospheric 
velocity fields as drivers for MHD waves. Three of the five drivers 
considered were spiral type drivers, based on observations of spiral 
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motions in the solar atmosphere (Bonet et al. 2008 [ 

Wedemeyer- 

Bohm & Rouppe van der Voort 2009, Bonet et al. 2010 

Wedemeyer 

et al.|2013|), these motions were modelled as circular, Archeme- 


dian and logarithmic motions. It was concluded that the logarithmic, 
Archemedian and uniform spiral drivers all generate similar (+10%) 
excited energy fluxes. The spiral expansion factors were selected 
arbitrarily in |Mumford et al. 1(2015) . This work analyses the effects 
of the spiral expansion factor on the MHD waves generated by the 
logarithmic spiral driver, motivated by the observational studies and 
constraints of |Bonet et al. 1(2008) . In |Bonet et al.| ([2008) magnetic 
bright points (MBPs) were observed spiralling in an inter-granular 
lane, where cold plasma sinks down into the convection zone. |Bonet| 
|et al.| ( [2008| ) fit the observed locations of the MBP with time to the 
equation for a logarithmic spiral, shown in Equation ([TJ, 

0 = -jr- In(r/fl), (1) 

where r is the radius of the spiral and a is a positive real constant, 
and obtained a value of B]~ l = 6.4 + 1.6 or B L = 0.15 for the 
dimensionless expansion factor parameter. 

In |Bonet et al.| ( |2010| > a larger sample of photospheric vortices 
were studied, despite not fitting spirals to the observed motions, 
a number density of photospheric vortices was calculated as d ^ 
3.1 X 10 -3 vortices Mm -2 minute -1 , which therefore provides an 
upper limit of the number of logarithmic spiral-like vortices in the 
solar photosphere. 

In this work we investigate the role of the spiral expansion 
factor ( B l ) in the generation of MHD waves in a non-potential 
Gaussian magnetic flux tube, embedded in a realistic stratified solar 
atmosphere. The observational result of |Bonet et al.| ( |2008) is used 
as a starting point and values ±3x and ±10x that value are then 
employed to give five points in the parameter space, centred around 
their result, which is illustrated in Figure [T] 


2 SIMULATION CONFIGURATION 


X X i- X - X X 

0.015 ' * 0.050 0.150 ' 0.450" " " " 1.500' 

B l 

Figure 1 . The parameter space of Bl used in this work shown as blue 
crosses, with the x-axis on a logarithmic scale. The green error bars show 
the fit uncertainty of the value observed by |Bonet et al.H2008) . 

photosphere. The domain is divided up into 128 3 grid points giving 
a physical size of 15.6 x 15.6 x 12.5 km 3 for each grid cell. All 
of the boundary conditions are open and therefore allow almost all 
non-linear perturbations to escape without significant reflection. 

The magnetohydro static background is perturbed during the 
simulations using a 3D Gaussian weighted logarithmic spiral ve¬ 
locity driver, as described by Equation (3]> ( [Mumford et al.|20l5) : 


Vx = A cos(e + d,) e -( 

Az 2 Ax 2 Ay 2 ) gin (2/r— ) , 

(3a) 

V* 2 + y 2 


A sin(e + </,) e 

( AA + Ax 2 + Ay 2 ) sin (27T-), 

(3b) 

V* 2 + y 2 

\ P) 

(3c) 


where: 

6 - tan~ l j, 0 = tan~ l > 

A = Ax = Ay = 0.1 Mm and A z = 0.05 and P = 180 s. Here, 
B l is the logarithmic spiral expansion factor discussed in Section]]] 
Figure [2] shows the calculated velocity profiles for the peak 
vertical height of the driver (z = 100 km). Overplotted on these 
profiles are streamlines that trace a logarithmic spiral with different 
expansion factors. 


The simulations performed for this study utilise a realistic stratified 
solar atmosphere constructed by taking the VALIIIc ( [Vemazza et al.| 
|1981| hydrodynamical properties and adding a non-potential self¬ 
similar magnetic field. The self-similar magnetic field configuration 
derived from the ones employed by Fedun et"ah] (2011| l and 


recently analytically described in|Gent et al. 


on|Schlut er & Te mesvary| ([T958); 


Deinzer 


2014k; based 
Low (1 


980k; 


Schtissler & Rempel ( 2005]), and identical to the one in Mumford 


et al. ( |2015] r A magnetic field is constructed via this method, then 


added to the hydrostatic background and then the pressure balance 
is satisfied using magneto-hydrostatic equilibrium as described by 
Equation ([2j, i.e. 


- (B b • V)B b + V ( 5^- | + V/j = pg. 


( 2 ) 


where B b is the background magnetic field, p is the density, and 
p is the pressure. Equation corrects the missing negative term 
in [Mumford et al.| ( |2015| ), the calculations are not affected as this 
was a typo. By using a magnetic footpoint strength of 120 mT and 
the background atmosphere as specified by the VALIIIc model, the 
resulting numerical domain has the plasma J3 > 1 at every point. 


3 ANALYSIS 

To quantify the MHD wave modes generated by the logarithmic 
spiral velocity drivers it is necessary to quantify the relative pro¬ 
portion of the excited MHD wave modes. The modes present in 
the domain are assumed to be uniquely determined by the three 
wave modes present in a uniform homogeneous plasma, namely, the 
fast magnetoacoustic mode, the slow magnetoacoustic mode and 
the Alfven mode. The numerical domain used in this work has the 
plasma J3 > 1 everywhere, therefore, we consider wave propagation 
in this regime. Under these conditions the three MHD wave modes 
are separable into three vector components of perturbation with 
respect to the magnetic field. The fast magnetoacoustic mode is the 
dominant mode in the parallel vector component with respect to 
the magnetic field. The slow magnetoacoustic mode is the domi¬ 
nant contributor to the vector component perpendicular to both the 
magnetic field vector and to the magnetic flux surfaces. The Alfven 
mode can be identified in the third vector component, found via 
the cross product of the parallel and perpendicular vector, parallel 
to the flux surface but perpendicular to the magnetic field vector. 
However, plasma geometry and conditions in the simulation domain 


The Sheffield Advanced Code (SAC) (Shelyag et al. 

2008 

i used 

make this approximation somewhat imperfect, because there are no 

in this work is configured identically to Mumford et al. ( 

2015 

). The 

clear MHD eigenmodes due to the physical coupling of the waves. 

domain has a spatial extent of 2.0 x 2.0 x 1.6 Mm 3 in x, 

y and 

Further, these three modes become degenerate in cylindrical geom- 


z respectively, with the origin in the z direction 61 km above the etry giving rise to sausage, kink, and fluting modes. Also, due to 
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(a) B l = 0.015 
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Figure 2. Cuts in the [ x-y ] plane through the driving velocity field. The normalised velocity is plotted in blue with velocity vectors overplotted in black and a 
streamline seeded at the centre plotted in red. A plot is shown for each value of Bi used in a simulation. 


the complex plasma conditions in the simulation domain the modes 
may become physically coupled meaning that it is impossible to 
completely separate the modes. Despite these complications the 
description of the modes based on the three vector components in 
the magnetic field frame is taken as a good way to describe, identify 
and quantify the MHD wave modes in the system. 

To identify theses waves via the vector components relative 
to the magnetic field the identification of a vector perpendicular to 
the magnetic field vector is required. In a 2D system this is a trivial 
step, however, in a 3D simulation it is ill-defined. The solution 
to this problem, used in this work, is to define a magnetic flux 
surface which encapsulates a constant amount of magnetic flux at 
all heights in the domain. This method is described in more detail 
and utilised in Mumfor d et al.| ( [201 5 1 >. The surface then allows 
the computation of a vector perpendicular to it and, thus, to the 
magnetic field lines it is constructed from. These ‘flux surfaces’ 
are initially constructed from a ring of axisymmetric field lines 
computed in the static background conditions. The field line seed 
points then move with the plasma velocity throughout the simulation, 
which results in the flux surface being constructed from the same 
field lines at all times in the simulation. The combination of the 
surface normal vectors and the magnetic field vector provide the 
information required to calculate the azimuthal vector via the cross 
product, which provides a third vector parallel to the surface but 
perpendicular to the magnetic field. 

These surfaces are constructed, using the VTK librar^ for 
three different characteristic initial radii (measured at the top of the 
domain) of 156 km, 468 km and 936 km from the centre of the 
domain, for each simulation, giving a good sampling through the 
differing plasma properties of the domain. This allows the analysis 

1 Visualistation ToolKit 5.10.0 (www. vtk. org) 


of the excited modes at different points in the domain, giving an 
overall picture of the waves. 

Using the flux surfaces, defined above, we can now decompose 
any vector quantity in the domain into the parallel, perpendicular 
and azimuthal components, allowing study of the velocity and mag¬ 
netic field perturbation vectors. While the velocity and magnetic 
field perturbation vectors are good for identifying and studying wave 
behaviour itself, to quantify the amount of each wave mode gen¬ 
erated the wave energy flux is computed using Equation 0 from 
[Bogdan et al.[ ( |2003| ). 

F wave = W + — (B* • B) V - — (v • 5) B^, (4) 

where subscript b represents a background variable, tilde represents 
a perturbation from the background conditions and p k represents 
kinetic pressure. 

The wave energy flux Equation 0 is decomposed onto the flux 
surface in the same way as the velocity vector, subject to the same 
limitations as the velocity. 


3.1 Results 

To assist in the visualisation and analysis of the results provided by 
the flux surfaces the vector components, for both velocity and wave 
flux, along one field line are extracted for all time steps and plotted 
as time-distance diagrams in Figures [3] and|4] 

Combining the decomposed velocity vector plotted in Figure 
[3] and the decomposed wave flux vector plotted in Figure]?] we can 
reliably describe the nature of the waves generated in the simulations. 
Overplotted on all panels in Figures [3] &[?] are the phase speeds for 
the background conditions, the dot-dashed line is the fast speed Vf, 
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the dashed line is the sound speed c s , the dotted line is the Alfven 
speed v A and the solid line is the slow speed v s . By comparing these 
characteristic wave mode speeds to the ridges in the time-distance 
diagrams it can be seen that in the panels for the torsional component 
(third panel in each figure), the dominant perturbation travels with 
the Alfven speed (solid line). We interpret this perturbation as an 
Alfven wave. For the perpendicular component (second panels) it 
can be seen that the dominant perturbation travels with the slow 
speed (solid line), therefore this perturbation is interpreted as a slow 
sausage mode. We can infer that this perturbation is likely to be a 
sausage mode perturbation due to the nature of the driver, in that it 
should not perturb the axis of the flux tube and, that we observe no 
significant displacement on the flux surfaces during the simulation. 
The most interesting result is shown for the parallel component (top 
panel in each figure), where for lower values of B L , the amplitudes 
are low, but the perturbations that are present travel with the slow 
speed (solid line). However, as B L increases the perturbations change 
form. There seems to appear a second, superimposed perturbation 
travelling with a speed close to that of the fast (or sound) speeds, 
which could be a fast sausage mode. This second perturbation seems 
to grow proportionally to B L , and can be seen to be dominant in 
Figure s[4d| and [4e| 

The wave flux graphs in Figure]?] are components normalised 
to the magnitude of the wave flux vector, thus showing the relative 
strengths of the components. Taking Figure [4a| for the B L = 0.015 
spiral it can be seen that most of the excited wave flux is in the 
azimuthal component, associated with the Alfven wave. As the 
expansion factor (. B L ) increases, the driver becomes more radial, and 
the flux starts to shift from the azimuthal component into the parallel 
component. This is interpreted as a change of the dominant mode 
from the torsional Alfven wave into a sausage mode with dominant 
velocity perturbations parallel to the field lines. Considering the 
range of B L , found by|Bonet et al.| ( |2008| » and illustrated in the range 
spanned by Figure |3c| and [3d] it can be seen that even within this 
parameter range the parallel component becomes substantially more 
dominant, meaning the change in spectrum of excited MHD wave 
modes is sensitive to the expansion factor of a spiral driver. 

|Mumford et al.| ( |20T5] > reported that, for the spiral drivers, there 
is a significant percentage of the wave energy flux contained in 
the perpendicular component. This appears to be inversely coupled 
to the spiral expansion factor of the driver, as it decreases propor¬ 
tionally with the azimuthal wave flux component. The size of the 
perpendicular component is also inversely proportional to the initial 
radius of the flux surface, as can be seen by its decrease in the three 
panels of Figure [5] 

This change in excitation of MHD waves is summarised in 

F 2 

1 || j_ a 

Figure 5 where the average value of —- ’-j- -j for all time is 

L ._. F a +F ± + F « 

plotted7"Figure[5] shows that, between the values of B L = 0.15 and 
B l = 0.45 there is a turning point where the torsional component 
becomes less dominant, with expansion factors larger than 0.15 
having the parallel component being the dominant component. This 
turning point occurs within the range of the fitted spirals in |Bonet| 
[et al.| ( |2008| > and, therefore, implies that photospheric spirals may 
indeed generate a variety of different MHD modes with varying 
strengths. 
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Figure 5. Comparison of percentage square wave energy flux for all simula¬ 
tions and all calculated flux surfaces. The parallel component of wave energy 
flux is shown as blue circles, the perpendicular component as green dashes 
and the azimuthal component as red crosses. The green shaded region shows 
the fit uncertainty in the value observed in |Bonet et al. 1(2008) 


MHD waves in a magnetic flux tube. The expansion factor of the 
logarithmic spirals was varied around a statisti cally determined ob- 

(2008). The 


Bonet et al. 


servational value of B = 6.4 + 1.6 from 
excited MHD waves were analysed using ‘flux surfaces’ computed 
from the magnetic field lines, from which physical vector parame¬ 
ters, such as velocity, were decomposed into a reference frame of 
parallel, perpendicular and azimuthal to the magnetic field lines. 
This decomposition enabled the identification of the excited MHD 
wave modes, and from computing the wave energy flux the relative 
strengths of the components was calculated. 

The average wave flux analysis for all time was used as an 
indication of the relative strength of the three components, in the 
magnetic field frame. The analysis indicates that between B L - 0.15 
and B l - 0.45 the torsional component, associated with the Alfven 
mode, becomes weaker than the parallel component, which becomes 
dominant. Assuming that the expansion rates of these observed 
spirals will be distributed over a range of possible values, they may 
not be generating the quantity of Alfven waves and fluxes previously 
assumed. 


4 CONCLUSIONS 
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Figure 3. Velocity time-distance diagrams for all simulated values of Bi for the surface with an initial top radius of 468 km. Shown in green are the phase 
speeds for the background conditions, the dot-dashed line is the fast speed Vf, the dashed line is the sound speed c s , the dotted line is the Alfven speed va and 
the solid line is the slow speed v s . Note that plasma ft > 1 for all heights in the domain. 
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Figure 4. Normalised wave energy flux time-distance diagrams for all simulated values of Bl for the surface with an initial top radius of 468 km. Shown in blue 
are the phase speeds for the background conditions, the dot-dashed line is the fast speed Vf, the dashed line is the sound speed c s , the dotted line is the Alfven 
speed va and the solid line is the slow speed v s . Note that plasma [3 > 1 for all heights in the domain. 
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